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ABSTRACT 



We report results from a Far Ultraviolet Spectroscopic Explorer (FUSE) survey of in- 
terstellar molecular hydrogen (H2) along 45 sight lines to AGN at high Galactic latitudes 
(|6| > 20°). Most (39 of 45) of the sight lines show detectable Galactic H2 absorption 
from Lyman and Werner bands between 1000 and 1126 A, with column densities ranging 
from Nh2 = 10 1 • 17_ 9,82 cm -2 . In the northern Galactic hemisphere, we identify many 
regions of low N H2 (< 10 15 cm" 2 ) between £ = 60° - 180° and at b > 54°. These "H 2 
holes" provide valuable, uncontaminated sight lines for extragalactic UV spectroscopy, 
and a few may be related to the "Northern Chimney" (low Na I absorption) and "Lock- 
man Hole" (low Nhi)- A comparison of high-latitude H2 with 139 OB-star sight lines 
surveyed in the Galactic disk suggests that high-latitude and disk H2 clouds may have 
different rates of heating, cooling, and UV excitation. For rotational states J = and 1, 
the mean excitation temperature at high latitude, (TqI) = 124±8 K, is somewhat higher 
than in the Galactic disk, (T$_ sk ) = 86 ±20 K. For J > 2, the mean (T cxc ) = 498 ±28 K, 
and the column-density ratios, N(3)/N(l), N(4)/N(0), and N(4)/N(2), indicate a com- 
parable degree of UV excitation in the disk and low halo for sight lines with Nh2 > 10 18 
cm -2 . The distribution of molecular fractions at high latitude shows a transition at 
lower total hydrogen column density (log N^ 1 ~ 20.38 ± 0.13) than in the Galactic disk 
(log Ng sk « 20.7). If the UV radiation fields are similar in disk and low halo, this 
suggests an enhanced H2 (dust-catalyzed) formation rate in higher-density, compressed 
clouds, which could be detectable as high-latitude, sheetlike infrared cirrus. 
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1. INTRODUCTION 

Molecular hydrogen (H2) is the most abundant molecule in the universe, comprising the ma- 
jority of mass in interstellar molecular clouds that eventually form stars. In the diffuse interstellar 
medium, with visual extinctions Ay < 1 mag, the absorbing clouds exhibit detectable H2 lines with 
molecular fractions, /h2, ranging from 10~ 6 at low H I column density up to ~40% (Savage et al. 
1977; Spitzer & Jenkins 1975; Shull et al. 2005). In sight lines with even greater extinction, the 
so-called translucent sight lines with Ay = 1 — 5 mag (van Dishoeck & Black 1986), the molecular 
fraction can be as high as 70% (Rachford et al. 2002). Even though H2 plays an important role 
in interstellar chemistry, many questions remain about its distribution, formation, and destruction 
in protostellar clouds. With the 1999 launch of the Far Ultraviolet Spectroscopic Explorer (FUSE) 
satellite, astronomers regained access to far ultraviolet (FUV) wavelengths below 1126 A, needed 
to study H2 via resonance absorption lines from its ground electronic state, X to the excited 
states, B 1 S+ (Lyman bands) and C 1 H U (Werner bands). In absorption by the cold clouds, the 
observed transitions originate in the ground [v = 0) vibrational state from a range of rotational 
states, primarily J = (para-H2 with total nuclear spin 5 = 0) and J = 1 (ortho-H2 with nuclear 
spin 5 = 1). Weaker absorption lines from excited states, J > 2, are also detected, which are 
produced by UV fluorescence and cascade following photoabsorption in the Lyman and Werner 
bands; see reviews by Spitzer & Jenkins (1975) and Shull & Beckwith (1982). 

From the FUSE mission inception, H2 studies have been an integral part of the science plan. 
The FUSE satellite, its mission, and its on-orbit performance are described in Moos et al. (2000) 
and Sahnow et al. (2000). The initial FUSE studies of H2 appeared among the Early Release 
Observations (Shull et al. 2000; Snow et al. 2000). More recent observations include FUSE surveys 
of H2 abundance in translucent clouds (Rachford et al. 2002) and diffuse clouds (Shull et al. 2005) 
throughout the Milky Way disk, in interstellar gas toward selected quasars (Sembach et al. 2001a, 
2004), in intermediate velocity clouds (IVCs) and high velocity clouds (HVCs) (Sembach et al. 
2001b; Richter et al. 2001, 2003); in the Monoceros Loop supernova remnant (Welsh, Rachford, &: 
Tumlinson 2002), and in gas throughout the Large and Small Magellanic Clouds (Tumlinson et al. 
2002). 

Sightlines toward active galactic nuclei (AGN) through the Galactic halo provide special probes 
of diffuse interstellar gas, including large-scale gaseous structures, IVCs, and HVCs. In contrast, 
surveys of OB stars in the Galactic disk typically have much higher column densities, Nh2 > 
10 18 cm -2 , except toward nearby stars too bright to be observable with the FUSE detectors. 
Therefore, the AGN halo survey is able to probe fundamental properties of H2 physics, tied to 
observable quantities such as metal abundances, H I, and IR cirrus. The halo survey also probes 
H2 in a new, optically thin regime (Nh2 < 10 17 cm~ 2 ) inaccessible to FUSE in the Galactic disk. 
These low-column clouds can provide important tests of models of H2 formation, destruction, and 
excitation in optically-thin environments where FUV radiation dominates the J-level populations. 
The targets in this survey were typically not observed for the primary purpose of detecting H2. 
Rather, these high-latitude AGN were observed to study Galactic O VI, for the Galactic D/H 
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survey, as probes of the interstellar medium (ISM) and intergalactic medium (IGM), and for their 
intrinsic interest in AGN studies. The relatively smooth AGN continua allow us to obtain clean 
measurements of interstellar absorption. 

Molecular hydrogen often contaminates the UV spectra of extragalactic sources, even when 
they are used for other purposes, such as studies of intervening interstellar and intergalactic gas. 
Because H2 has such a rich band structure between 912-1126 A, one must remove H2 Lyman and 
Werner band absorption by using a model for the column densities, N(J), in H2 rotational states 
J = — 5 and occasionally even higher. This H2 spectral contamination becomes detectable at 
log Nh2 > 14. It becomes serious at log Nh2 > 16, when the resonance lines from J = and J = 1 
become strong and selected lines from higher rotational states blend with lines from the ISM and 
IGM. A few H2 absorption lines blend with key diagnostic lines (O VI AA1031.926, 1037.617) and 
occasionally with redshifted IGM absorbers. Our AGN sample is useful for another reason: we 
are able to assess the range of H2 column densities and corrections to FUV spectra taken to study 
the ISM and IGM. This subtraction was performed in studies of intergalactic O VI (Savage et al. 
2002; Danforth & Shull 2005), of high-latitude Galactic O VI (Wakker et al. 2003; Sembach et al. 
2003; Collins, Shull, & Giroux 2004), and of D/H in high velocity cloud Complex C (Sembach et al. 
2004). 

Observations of the hydrogen molecule, its abundance fraction in diffuse clouds, and its rota- 
tional excitation provide important physical diagnostics of diffuse interstellar gas (Shull & Beckwith 
1982). For example, the rotational temperature, Tqi, of the lowest two rotational states, J = and 
J = 1, is an approximate measure of the gas kinetic temperature in cases where collisions with H + 
or H° are sufficiently rapid to mix the two states. In diffuse clouds, the H2 rotational temperature, 
(Toi) = 77 ± 17 K, measured by Copernicus (Savage et al. 1977) provides one of the best indicators 
of the heating and cooling processes in the diffuse ISM within 1 kpc of the Sun. Our FUSE survey 
(Shull et al. 2005) extends this study to more distant OB stars in the Galactic disk and finds similar 
temperatures (Tqi) = 86 ± 20 K. 

For higher rotational states, the level populations are characterized by an excitation temper- 
ature, T exc , that reflects the FUV radiation field between 912-1126 A, which excites and pho- 
todissociates H2. Populations of the J > 2 rotational levels in the ground vibrational and ground 
electronic state are governed by absorption in the H2 Lyman and Werner bands. These electronic 
excitations are followed by UV fluorescence to the ground electronic state and infrared radiative cas- 
cade through excited vibrational and rotational states, often interrupted by rotational de-excitation 
by H°-H2 collisions. In approximately 11% of the line photoabsorptions, for a typical ultraviolet 
continuum, H2 decays to the vibrational continuum and is dissociated. This peculiarity of H2 al- 
lows the molecule to "self-shield" against dissociating FUV radiation, as the Lyman and Werner 
absorption lines become optically thick. Detailed models of these processes (Black k, Dalgarno 
1973, 1976; Browning, Tumlinson, & Shull 2003) can be used to estimate the FUV radiation field 
and gas density. 
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In this paper, we present a survey of H2 toward 45 AGN at Galactic latitudes |6| > 20° (see 
Figure 1). To produce detectable H2, with /h2 > 10 -5 , these clouds typically need hydrogen gas 
densities uh > 10 cm -3 , so that the H2 formation rates are sufficient to offset UV photodissociation. 
At column densities Nh2 > 10 15 cm~ 2 , the molecules begin to self-shield from UV photodestruction. 
The selection of high-latitude AGN background targets allows us to minimize the amount of H 2 from 
the Galactic disk, although we cannot exclude it altogether. Our high-latitude survey accompanies 
a larger survey (Shull et al. 2005) of H2 toward 139 OB-type stars in the Galactic disk. These 
OB stars at < 10° typically show larger H2 column densities and greater molecular abundances 
than the AGN sight lines. Consequently, the AGN sight lines provide a better sample of diffuse 
H2-bearing clouds above the disk, and perhaps into the low halo. Various properties of the gas, 
including molecular fraction, rotational temperature, and connections with infrared cirrus clouds 
(Gillmon & Shull 2005), suggest that the high-latitude sight lines sample a different population of 
diffuse absorbers than the disk survey. The halo H2 survey may also be representative of sight lines 
that pass outside the inner disks of external galaxies. 

The organization of this paper begins with a description of the AGN sample and our data 
acquisition, reduction, and analysis (§ 2). This is followed in § 3 by discussion of the results on 
H2 abundances, molecular fractions, and rotational excitation. We examine both the rotational 
temperature, Tqi, for J = and J = 1, and the excitation temperature, T exc , for higher rotational 
levels. We also show the distribution of H2 in column density and the spatial distribution throughout 
the northern Galactic hemisphere, particularly the conspicuous "H2 holes" with Nh2 < 10 15 cm~ 2 
at b > 54°. We conclude in § 4 with a comparison of these high-latitude results to those from 
Galactic disk surveys by Copernicus (Savage et al. 1977) and FUSE (Shull et al. 2005). 

2. DATA ACQUISITION AND REDUCTION 

2.1. FUSE Observations 

Table 1 lists our target names and their Galactic coordinates, source type, B magnitude, color 
excess E(B-V), and observational parameters. We chose our targets based on their signal-to-noise 
ratio (S/N), source type, and Galactic latitude. All 45 of our targets are at \b\ > 20°, a subset of the 
219 targets selected in Wakker et al. (2003) as candidates for the analysis of Galactic O VI. Of the 
219 targets, Wakker et al. found 102 appropriate for O VI analysis, for which they binned the data 
by 5-15 pixels (the FUSE resolution element is 8-10 pixels or about 20 km s -1 ) in order to meet 
an imposed requirement of (S/N)bi n > 4. We set a similar S/N requirement for our H2 analysis. 
However, we exclude those targets that required binning by more than 8 pixels, because the loss 
of resolution caused by overbinning becomes problematic for the analysis of relatively narrow H2 
lines, in comparison with the typically broad (> 50 km s _1 ) O VI profiles. Of the 219 targets, 
55 met our requirement of (S/N)bi n > 4 or (S/N) p i x > 2 with 4-pixel binning. We excluded halo 
stars and restrict our sample targets to AGN plus a few starburst galaxies in order to minimize 
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the difficulties presented by stellar continua. We checked all targets that were excluded for low 
S/N to see whether additional FUSE data had been acquired since the Wakker et al. (2003) survey. 
The additional data increased the S /N for seven of these targets sufficiently to include them in our 
sample. 

All observations were obtained in time-tag (TTAG) mode, using the 30" x 30" LWRS aperture. 
The resolution of FUSE varies from R = A/AA = 15,000 - 20,000 across the far-UV band, and 
it can also vary between observations. However, our results rely on measured equivalent widths 
that are effectively independent of instrumental resolution. Data were retrieved from the archive 
and reduced locally using CalFUSE v2.4. The CalFUSE processing is described in the FUSE 
Observers Guide (http://fuse.pha.jhu.edu). Raw exposures within a single FUSE observation were 
coadded by channel mid-way through the pipeline. This can produce a significant improvement 
in data quality for faint sources (such as AGN) since the combined pixel file has higher S/N than 
the individual exposures, and consequently the extraction apertures are more likely to be placed 
correctly. Combining exposures also speeds up the data reduction time substantially. 

Because FUSE has no internal wavelength calibration source, all data are calibrated using 
a wavelength solution derived from in-orbit observations of sources with well-studied interstellar 
components (Sahnow et al. 2000). This process leads to relative wavelength errors of up to 20 
km s _1 across the band, in addition to a wavelength zero point that varies between observations. 
For targets with data from multiple observations, the spectra from separate observations were 
coaligned and coadded to generate a final spectrum for data channels LiFla and LiF2a (the other 
six data channels were not used in our analysis). The coalignment was performed on one single line 
per channel: Ar I at 1048.22 A for LiFla, and Lyman (1-0) P(3) at 1099.79 A for LiF2a. Since we 
did not correct for possible relative wavelength errors across the wavelength range before coadding, 
our method of using a single line to perform the alignment may blur out weak lines at wavelengths 
well away from 1048.22 A (LiFla) and 1099.79 A (LiF2a). 

Once the final coadded spectra have been obtained, there may still be errors in the wavelength 
solution. Because the 21 cm spectra are acquired with ground-based telescopes, their wavelength 
solution is more accurate than that of FUSE data. We attempted to correct the FUSE wavelength 
solution by plotting the 21 cm emission line and Ar I 1048.22 A absorption line in velocity space 
and shifting the FUSE spectra to align the centroids. Although it is not necessarily the case that 
the neutral gas (H I) must align in velocity with low- ionization metals (Ar I, Si II, Fe II) or with H2, 
it is a plausible assumption that provides a standard technique for assigning the FUSE wavelength 
scale. Although relative wavelength errors may exist across the full wavelength range, we place 
our trust in the wavelength solution near 1048 A and use H2 lines in the (4-0) and (5-0) Lyman 
bands between 1036-1062 A to determine the velocity of the H2 clouds (see § 2.2). The S/N of the 
co-added data ranges from 2-11 per pixel, and it also varies with spectral resolution, which is not 
fixed in our survey. Most of the data were binned by 4 pixels before analysis, with the rare case of 
binning by 2 or 8 pixels. 
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2.2. Data Analysis 

In our analysis of H2, we use the signature absorption bands in the FUSE LiF channels between 
1000 and 1126 A. These bands arise from the Lyman and Werner electronic transitions from the 
ground (v = 0) vibrational state and rotational states J = — 7. Wavelengths, oscillator strengths 
(/), and damping constants for these lines were taken from Abgrall et al. ( 1993a, b). Because of 
the rapid decrease in FUSE sensitivity and decreased S/N in the SiC channels at wavelengths 
less than 1000 A, we restrict our analysis to ten vibrational-rotational bands, Lyman (0-0) to 
(8-0) and Werner (0-0). Additional information can be gleaned from higher vibrational bands 
(A < 1000 A) in the SiC channels, but we do not present this data here. To maximize efficiency, 
we restricted our analysis to the LiFla and LiF2a channels, which provide the best combination of 
wavelength coverage and sensitivity. These channels provide higher S/N than the SiC channels in 
this wavelength range and include the Lyman (2-0) and (1-0) bands. 

Seven AGN sight lines show H2 absorption in J = 4, and absorption lines from J = 5 were 
detected toward one target (HS 0624+6907). Approximately one-third (15) of the sight lines show 
substantial H2 column densities, ranging from log Njj2 = 18.0 - 19.8, with damping wings in the 
J = and J = 1 lines. Half the AGN spectra show weaker H2 absorption with log Nh2 ~ 14 - 17. 
A summary of our results appears in Table 2, where we list H2 column densities in each rotational 
state ( J = — 5) and the doppler parameter (b) from the CoG. 

Our analysis routines are adapted from Tumlinson et al. (2002), who developed our standard 
software for analyzing H2. These routines allow us to determine equivalent widths for as many 
H2 lines as possible, and then to apply a curve of growth (CoG) to arrive at column densities in 
individual rotational states, N(J). From this information, we can infer various properties of the 
gas. We fitted all detected H2 lines between 1000 and 1126 A that are not visibly compromised by 
continuum problems or line blending, using either a Gaussian or Voigt profile in order to measure 
the equivalent width, W\. For most sight lines, the fits to the absorption lines are carefully tailored 
on a line-by-line basis to choose the correct continuum and adjust for any line blending. When a 
proper fit to a line or group of lines has been obtained, the equivalent width W\, uncertainty in 
W\, and central wavelength are stored in an electronic table for use by the CoG fitting routines. In 
select cases, we fitted the J > 2 lines with a more automated routine, also described in Tumlinson 
et al. (2002), which queries the user at each wavelength corresponding to an H2 line and allows it 
to be fitted with a single Gaussian component. This routine is useful for targets with relatively flat 
continua and minimal line blending. The automated routine also stores W\, uncertainty in W\, and 
central wavelength in a table. The J = and 1 lines are fitted by hand with a Voigt profile since 
they may be strong enough to have non-Gaussian profiles. For lines between 1036-1062 A in the 
Lyman (5-0) and (4-0) bands, with well-determined centroids, we converted the central wavelength 
to a velocity using the Doppler formula and then averaged the velocities. We list these H2 cloud 
velocities (LSR) in Table 3, together with velocities of the various 21-cm emission components 
(Wakker et al. 2003). 
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A single sight line is likely to intercept many clouds containing H2, that may or may not be 
offset in velocity. One of the potential sources of systematic uncertainty is our inability to resolve 
components of H2 absorption by different clouds along the line of sight. At the FUSE resolution of 
20 km s _1 , we typically cannot separate absorption components in the "LSR core", from —20 to 
+20 km s _1 . Components separated by > 30 km s _1 can sometimes be resolved, depending on the 
line strengths and widths. We have 24 targets in common with Richter et al. (2003), who focused 
on H2 in intermediate velocity clouds with \v\ = 30 — 90 km s -1 . They report the existence of an 
IVC in 7 of our sight lines and cite a possible IVC in 6 others. Because we focus on low-velocity 
H2, we generally avoid IVCs in our line fitting. We have determined the column density for an IVC 
in three cases: 3C 273, Mrk 876, and NGC 4151. For 3C 273, we measured H 2 at +25 km s -1 , 
which we do not consider a true IVC. For Mrk 876, we fitted the IVC lines along with the LSR 
core lines (double component fits). For NGC 4151, the H2 lines arise mainly from IVC gas, but 
they may have contributions from an LSR component. Out of the 7 targets with confirmed IVCs, 
it is possible in one case that the IVC contributed to our measured LSR H2 lines. For the 6 targets 
with possible IVCs, it is possible in four cases that an IVC contributed to our fitted LSR H2 lines. 
We also note 4 targets where Richter et al. (2003) claim "no evidence" for an IVC, but we believe 
an IVC contributed to our measured H2 lines. Our comments on velocity components and IVCs 
are summarized in Table 3. Multiple unresolved components in the LSR core, as well as possible 
contributions from IVCs, may produce a composite CoG [b and N(J)} that reflects the properties 
of the different velocity components. 

Once the line fitting is complete, the software described in Tumlinson et al. (2002) produces 
an automated least-squares fit to a CoG with a single doppler b parameter for all J levels. For the 
two sight lines noted above (3C 273 and Mrk 876), we make a two-component fit for the LSR core 
and IVC. Table 2 gives the best-fit values of N(J) and single-6, with 1 a confidence intervals. For 
low-S/N sight lines or those with low column densities, it is sometimes difficult to measure enough 
lines to constrain the CoG. This results in large error bars on N(J) and b. In some cases, model 
spectra of varying N(J) and b are overlaid on the spectrum to narrow the range of reasonable 
values. In many cases, the Lyman (1-0) and (0-0) J = 3 lines are crucial in constraining the b 
value. Large error bars on N(J) and b often result when these lines are absent. 

For sight lines with no discernable H2 lines, our methodology is as follows. We place upper 
limits on the equivalent widths of the strong Lyman (7-0) R(0) and R(l) lines and convert these 
limits to column density limits assuming a linear CoG. For the (4 a) limiting equivalent width of 
an unresolved line at wavelength Ao, we use the following expression: 

w Yim = 1^0 ^ 

A (A/AA)(S/N) 1 ' 

Since the resolution of FUSE varies across the band, we conservatively assume R = A/AA « 15, 000 
for all upper limits. The limiting equivalent widths range from 8 mA for high-S/N data to 43 mA 
for low-S/N data. This corresponds to a limiting column density range of log Nh2 > 13.8 — 14.6. 
The upper limits are included in Table 2. 
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2.3. H I COLUMN DENSITIES 

In order to fully interpret our observations of interstellar H2, we obtain supplementary data 
on the atomic hydrogen (H I) column densities associated with the observed H2. Since our sample 
of 45 sight lines is a subset of that in Wakker et al. (2003), we use their compilation of the 21-cm 
spectrum, as well as the component structure and N(H I) determined by fits to the data. They 
collected 21-cm spectra from the following sources: Leiden-Dwingeloo Survey (35' beam), Villa 
Elisa telescope (34' beam), Green Bank 140-ft telescope (1' beam), and the Effelsberg telescope 
(9/7 beam). Wakker et al. (2003) provide the 21-cm spectrum, as well as the component structure 
and N(H I) determined by fits to the data. 

To derive Nhi, we must determine which of the H I components are physically associated with 
the observed H2. This assignment can be uncertain and somewhat subjective, so we take care to 
choose a "best value" and a range of H I column densities, from N min to N max . Figure 2 illustrates 
our technique for Mrk 509, with the 21-cm emission line, H2 Lyman (4-0) R(0) absorption line, and 
Ar I 1048.22 A absorption line aligned in velocity space. We shift the FUSE data with respect to 
the 21-cm data in order to align the H I emission and Ar I absorption as described in § 2.1. We 
then select as Nb es t the sum of column densities in the H I components that are associated with H2 
absorption, based on the coincidence of the radial velocities. Of these selected components, at the 
20 km s -1 resolution of FUSE, there remains some uncertainty about which are actually associated 
with H2. In order to reflect this systematic uncertainty and to be conservative with errors, we 
adopt a range for the H I column density, including a floor of systematic uncertainty of ±0.03 in 
log Nhi- For the lower limit, N m i n , we use the smallest individual column density that could be 
responsible for the H2 absorption within the velocity range of the R(0) or R(l) line widths. For the 
upper limit, N max , we use the sum of all the H I column densities that can associated kinematically 
with the H2. For sight lines with no discernable H2 lines, we adopt a lower limit on Nhi as the 
smallest of the components that is not classified as an IVC or HVC and that satisfies the condition 
Nhi > 10 19 cm' 2 . 

Further systematic uncertainty arises from the larger beam size of the 21-cm data compared 
with the pencil-beam sight lines of the H2 data. It has been shown that a smaller beam gives a better 
approximation to the H I column density in the pencil beam (Wakker & Schwarz 1991; Wakker 
et al. 2001, 2002). On the other hand, Lockman & Condon (2005), using the 9'.8 Green Bank 
telescope, find that Nhi is highly correlated spatially. In Table 4, we provide various information 
on H2 and H I, as well as the 21-cm beam size. 
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3. RESULTS 

3.1. General Results 

As noted earlier, H2 absorption is present in most of the high-latitude sight lines. Figure 3 
shows a sample FUSE spectrum of PG 1211+143, an AGN sight line that has been studied ex- 
tensively with both FUSE and Hubble Space Telescope (HST) by Penton, Stocke, & Shull (2004) 
and Tumlinson et al. (2005). This sight line has a substantial column density, log Nh2= 18-381^4, 
producing a clear H 2 band structure throughout the FUV. The lines are saturated, but strong 
damping wings are not yet present. Figure 4 shows a sequence of four FUSE spectra of AGN, in 
order of increasing H2 column density. These range from Ton S180, where we detect no H2 to a 
(4cr) limit log N H2 < 14.37 (< 13.98 in J = and < 14.14 in J = 1), to ESO 141-G55, in which 
log Nh2 = 19.32 ± 0.07 with prominent damping wings in the R(0) and R(l) lines arising from 
J = and J = 1. This montage shows a number of H2 absorption lines from J = — 4 in the 
important Lyman (4-0) band, along with interstellar lines of Ar I and Fe II, the former of which is 
used to define the LSR velocity scale. 

Table 2 lists the H 2 column densities in individual rotational states, J = — 5 and the doppler 
(6) parameters from the CoGs. Figure 5 displays the distribution of H2 column densities, from the 
detectable limit, log Njj2 ~ 14.2, up to the maximum observed value log Nh2 ~ 19.8. Except for 
the dip at log Nh2 = 17-18, the distribution is fairly flat in column density. Small-number statistics 
preclude any firm conclusion about whether this dip is real, but it does suggest two populations of H2 
absorbers. For each sight line, Table 4 lists the total H2 and H I column densities, molecular fraction 
(/H2), rotational temperature, T01, of J = and J = 1 states, and the excitation temperature, 
T cxc , that characterizes the higher-J states. These temperatures were derived by least-squares fits 
of the column densities, N(J), to the form, 

- gjexp -— . (2) 



Z H2 V kT y 

Here, gj = gs (2 J + 1) is the statistical weight of rotational level J, with spin factor g$ = 1 or 3 
for para- or ortho-H2, respectively, Ej is the excitation energy of level J, Zyo. is the H2 rotational 
partition function, and T denotes Tqi for rotational states J = and 1 or T exc for excited states 
J > 2. 



3.2. Molecular Fraction 

Among the parameters listed in Table 4 is the average molecular fraction, 

2N(H 2 ) 

Jm N(H I) +2N(H 2 ) ' { ' 

which expresses the fraction of all hydrogen atoms bound into H2 molecules. We define the total 
column density of hydrogen as Nh = Nhi + 2Nh2, derived from Nh2, the H2 column density in 
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all rotational states, and Nhi, the neutral hydrogen column density derived from 21-cm emission 
(§ 2.3). Error bars on log Nh2 are found by propagating the uncertainties on log N(0) and log N(l), 
which dominate the H2 column density. 

In optically thin clouds, the density of molecules can be approximated by the equilibrium 
between formation and destruction, 



In this formula, the numerical value for /h2 is scaled to fiducial values of hydrogen density, nn (30 
cm -3 ), H2 formation rate coefficient, R (10~ 17 cm 3 s -1 ), and mean H2 pumping rate in the FUV 
Lyman and Werner bands, [3q = 5 x 10~ 10 s _1 . The H2 photodissociation rate is written as (/diss)/?, 
where the coefficient (/diss) ~ 0.11 is the average fraction of FUV excitations of H2 that result in 
decays to the dissociating continuum. The H2 formation rate per unit volume is written as nunuiR, 
where R depends on the gas temperature, grain surface temperature, and gas metallicity (Z). The 
metallicity dependence comes from the assumed scaling of grain-surface catalysis sites with the 
grain/gas ratio. For sight lines in the local Galactic disk, this rate coefficient has been estimated 
(Jura 1974) to range from R = (1 — 3) x 10~ 17 cm 3 s -1 for solar metallicity. This standard value 
for R is expected to apply when H I can stick to grain surfaces at suitably low temperatures of gas 
(Tgas < 300 K) and grains (T gr < 100 K) as discussed by Shull h, Beckwith (1982) and Hollenbach 
& McKee (1979). Although some infalling halo gas may have metallicities as low as 10-25% solar, 
as in HVC Complex C (Collins et al. 2003), we expect that most of the H2-bearing clouds in the 
low halo will have near-solar abundances and grain/gas ratios. 

Once sufficient column density of H2 builds up, with Nh2 > 10 15 cm -2 , the cloud becomes 
optically thick in individual Lyman and Werner lines, and the dissociation rate begins to drop, 
owing to self-shielding in the lines. The molecular fraction rises and makes an abrupt transition 
to much higher values, /h2 > 10 -2 , when the strong R(0) and R(l) lines develop damping wings 
at Nh2 > 10 18 cm -2 . For our halo survey, Figure 6 shows the molecular fractions, /h2, vs. total 
H column density. The molecular fraction rises from low values, /h2 ~ l0 -5 5=t0 ' 5 , to values above 
10~ 2 , typical of sight lines in the disk. The transition for halo sight lines occurs over a range log 
Nh ~ 20.2 — 20.5 (we quote a transition at log Nh = 20.38 ± 0.13), approximately a factor of two 
below the similar transition in the Milky Way disk at log Nh = 20.7 (Savage et al. 1977; Shull et al. 
2005). The precise location of the halo transition is not well determined from just 45 sight lines. 
Some broadening of the distribution is expected, since these sight lines may sample a mixture of 
gas in the low halo (low /H2) and denser clouds in the disk (higher /h2)- 

The large observed fluctuations in H2 column density require a very patchy ISM in the halo. 
For 25 sight lines in the northern Galactic hemisphere at latitude |6| > 20°, the mean total H2 
column density is (Nn2sin(6)) = 16.80, with a range from 10 14 to 10 19 ' 8 cm -2 . In the Copernicus 
H2 survey, Savage et al. (1977) found an average mid-plane density no(H2) = 0.036 cm -3 for 76 
stars, which they increased to 0.143 cm" 3 after correction for sampling biases due to reddening. 
If the H2 absorption was distributed smoothly in an exponential layer, n(z) = no exp(— z//i), 
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with scale height h, the integrated column density toward a quasar at latitude b would be Nh2 
= noh/ sin(6). For h ~ 100 pc, the column density would exceed Nh2 = 10 19 cm" 2 , even at high 
latitude. Obviously, a smooth distribution of H2 is in disagreement with our observations. 

The obvious explanation for the wide variations in Njj2 is a clumpy medium. However, the 
observed range (and fluctuations) in Nh2 constrain the geometry and hydrogen density of the H2- 
bearing "clouds". For example, consider an ensemble of spherical clouds of radius r c \ = (3 pc)r3 
and density nu = (30 cm~ 3 )n3o. If H2 is in formation-destruction equilibrium with the mean FUV 
radiation field, the molecular density inside the cloud is: 

nm * KUjj * (L64 X 10 ~ 4 cm ~ 3 ) n ^-i7 (j) , (5) 

where, as before, we have scaled the H2 formation rate coefficient R to 10 -17 cm 3 s _1 . The typical 
H2 column density through one such cloud would be 

AT(H 2 ) = „ (2 x 1( p cm- 2 )n 2 ^ 17 (|) r 3 . (6) 

In order to explain the total H 2 column density, a sight line to an AGN would have to intercept 
hundreds of such absorbers. The fluctuations in the number of clouds along different sight lines 
would then be quite small, certainly much less than the range observed. On the other hand, a 
smaller number of dense, sheetlike clouds, with a patchy distribution like the infrared cirrus, would 
be consistent with the data. These clouds would each have more Nh2, so the required number 
of interceptions would be smaller. Some high-latitude sight lines would have very large column 
densities, while others would be quite low, as in the "H 2 Holes" that we discuss in § 3.5. We return 
to the cirrus model in § 3.4. 



3.3. Rotational Temperature 

Figure 7 shows the J = 0-1 rotational excitation temperature, T01, which can measure the 
kinetic temperature in dense clouds, where H°-H 2 and H+-H2 collisions dominate the mixing of 
the J = and J = 1 rotational levels. The solid squares indicate high-Nn2 sight lines, where 
models suggest that the rotational temperature should be close to the gas kinetic temperature. 
The quantity Tqi is obtained from the expression 

AE 01 /k 

01 H( 9l /g )N(0)/N(l)} ' 1 > 

where gi/go = 9 is the ratio of the statistical weights of the J = 1 and J = rotational levels, and 
AE i/k = 170.5 K is the excitation temperature of the J = 1 level. 



To assess the statistical properties of the H2 rotational temperature in the halo, we selected 29 
sight lines with sufficiently small error bars on N(0) and N(l) to provide useful information on Tqi. 
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The halo rotational temperatures range from Tqi = 68-252 K, with a median of 139 K, a range 
that extends higher than the FUSE Galactic disk values, all but two of which lie between 55-120 K 
(Shull et al. 2005). For these 29 sight lines, the mean rotational temperature, {Tqi) = 124 ±8 K, 
is larger than the values found in the Copernicus H2 survey (77 ± 17 K, Savage et al. 1977) and in 
our FUSE Galactic disk survey (86 ± 20 K, Shull et al. 2005). For the 15 sight lines with log N H2 
> 17, we find (T 01 ) = 109 ± 7 K. 

An elevated kinetic temperature in the low halo could signal a shift in thermal balance between 
UV photoelectric heating by dust grains and radiative cooling, primarily by the [C II] 158 /im fine 
structure line. Both the heating and cooling rates should depend on the gas-phase metallicity, if 
the grain abundance scales with metallicity. In the low halo, we expect 0.1 < (Z/Z & ) < 1.0, with 
most clouds showing 30%-100% solar metallicities. Even though the abundances for these clouds 
might differ, the equilibrium gas temperature should be nearly independent of metallicity (Shull & 
Woods 1985), owing to cancelling dependences of heating and cooling on Z. The distribution of 
measured Tqi is broad, and some of the sight lines have large error bars. The median temperature 
is 121 K, although a minimum-x 2 fit gives 83 ± 6 K. Therefore, it may be premature to quote a 
higher temperature of these clouds. Additional studies of high-latitude sight lines would be most 
helpful. 



3.4. Excitation Temperature 

Figure 8 shows the higher- J excitation temperature, T exc , which reflects non-thermal UV 
fluorescent pumping of the high-J levels. The mean excitation temperature in the Galactic halo 
survey is (T cxc ) = 505 ± 28 K, with a range from 200-852 K. The median value is 501 K, while a 
X 2 fit to a constant temperature gives 416 ± 15 K. These values are similar to those seen by the 
Copernicus survey (Spitzer & Jenkins 1975; Shull & Bcckwith 1982), but slightly higher than that 
(326 ± 125 K) found in our FUSE disk survey (Shull et al. 2005). 

Figure 9 shows three ratios of high- J rotational states: N(3)/N(l), N(4)/N(2), and N(4)/N(0) 
versus log Nh2 in both low-latitude and high-latitude H2 surveys. The high-latitude survey contains 
32 sight lines with J = 3 measurements. Of these, only 8 sight lines have measured values of both 
N(3) and N(4), and all have substantial column densities, log Nh2 > 18. The J-ratios in these 
absorbers reflect the relative effects of UV excitation and collisional de-excitation for para-H2 
( J = and 2) and ortho-H2 ( J = 1 and 3), respectively. For the N(4)/N(2) ratio, the approximate 
ranges for the surveys are: 

• High-Latitude: log N(4)/N(2) = -2.8 to -1.8 and log N H2 = 18.5-20.0 

• Low-Latitude: log N(4)/N(2) = -4.0 to -1.0 and log N H2 = 19.5-21.0 

In general, these two distributions appear to cover similar ranges, particularly for absorbers 
with log Nh2 > 18. The greatest distinction between these populations occurs for the ratio 
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N(3)/N(l) in a subset of 16 high-latitude sight lines with low H2 column density, log Nh2 = 14-17, 
and have ratios log N(3)/N(l) = —1.5 to —0.5, well above the mean value in the disk (—3.0). These 
clouds are more optically thin in the dissociating FUV radiation, and thus have not built up large 
column densities in J = or 1. For the 8 sight lines with data for both J = 3 and J = 4, the degree 
of UV excitation appears similar in both latitude regimes. This result is consistent with models of 
the Galactic radiation field above the disk (Wolfire et al. 1995; Bland-Hawthorn & Maloney 1999), 
which suggest that the FUV radiation field drops by only ~ 50% at elevation z « 1 kpc. 

The high-latitude sight lines typically have a factor of 30 lower H2 column densities than those 
in the disk. For the disk absorbers, the highly excited J populations probably reside in surface 
layers surrounding a core with greater column densities in N(0) and N(l). Our recent models of 
H2 formation and excitation processes (Browning et al. 2003) can be used to estimate the effects 
UV radiation field and gas density from ratios such as N(4)/N(2) or N(3)/N(l). Figure 10 shows 
models that confirm these predictions in our high-latitude data. 

Taken together, these ratio distributions suggest that the degree of rotational excitation in the 
high-latitude sight lines is comparable to that in our disk survey (Shull et al. 2005) for high-Nn2 
sight lines. At first glance, this conclusion may seem inconsistent with the observed distribution of 
molecular fractions, fm, with column density, Nr. For the halo sight lines, the transition in /h2 
occurs at 50% lower Nh- As shown in equation (4), the molecular fraction scales as /h2 °c (n#i?//3). 
Thus, a shift in the distribution could arise either from a reduced rate of H2 dissociation (low (3) or 
an increased rate of H2 formation, nxR, on grain surfaces. The formation rate could be increased 
by a higher grain abundance or surface area (unlikely at the metallicities of halo clouds) or by 
a higher density nu resulting from cloud compression. Because we cannot supply clear evidence 
for an anomalously efficient population of dust grains in these high-latitude clouds, we favor the 
high-n# explanation. However, additional studies of the dust properties of these absorbers would 
be useful. 

Therefore, our best explanation of the observed J-ratios and the shift in the /b.2(^h) transition 
is that comparable UV dissociation rates in the halo and disk are offset by enhanced grain formation 
of H2 arising from halo clouds with higher density than in the disk. For example, these halo H2 
absorbers may be compressed sheets associated with the IR cirrus, as suggested for the strong H2 
toward ESO 141-G55 (Shull et al. 2000). We have followed up this idea (Gillmon k Shull 2005), 
where we find a clear association between H2 and IR cirrus clouds. Specifically, we see the same 
self-shielding transition of /h2 versus 100 fj,m cirrus intensity as seen with Nh- This suggests that 
the cirrus clouds with Nh2> 10 18 ' 5 cm" 2 have a substantial molecular content, ranging from 1-30%. 

3.5. Holes in the H2 Distribution 

Many of our high-latitude AGN targets lie behind regions of low H2 column density. As noted 
earlier, most of our sight lines show Galactic H2 absorption with detected column densities ranging 
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from Nh2 = 10 14,17-19 - 67 cm -2 . In the northern Galactic hemisphere, 8 sight lines have weak H2 
absorption (N H 2 < 10 15 cm -2 ) located between / = 60° - 180° at b > 54°. We have identified 
these high-latitude regions as "H2 holes", a patchy network of interstellar gas with both low 21-cm 
emission and low H2 absorption (Figure 11). 

The 8 sight lines through the H2 holes provide special locations to probe diffuse halo gas, 
including large-scale halo structures, IVCs, and HVCs. The structure of these holes may have 
been shaped by infalling H I clouds, as well as by outflowing gas from the "Northern Chimney", 
an interstellar cavity inferred from low Na I absorption at z < 300 pc (Lallement et al. 2003). 
Large portions of the holes have Nh2 < 10 15 cm~ 2 . These regions of low Nh2 are analogous to the 
"Lockman Hole" of low Nhi (Lockman et al. 1986), but they have much greater extent. Figure 
11 compares our H2 survey with the distribution of Nhi in the northern Galactic hemisphere at 
b > 45°. 

4. DISCUSSION AND SUMMARY 

Our FUSE high-latitude survey was able to probe molecular fractions below and above the 
transition value at log Nh ~ 20.38 ±0.13. We compared this transition to that seen in the Galactic 
disk surveys (Savage et al. 1977; Shull et al. 2005). Our sample, and possibly an enlarged future 
survey, can be used to examine fundamental properties of H2, tied to O VI, H I, and IR cirrus, 
in an optically thin regime (log Nh2 < 17) inaccessible to FUSE in the Galactic disk. These low- 
column clouds provide important tests of our theoretical models of H2 formation, destruction, and 
excitation in optically-thin environments, where FUV radiation dominates the J-level populations. 

These high-latitude sight lines are also relevant for far-UV studies of interstellar gas, inter- 
galactic gas, and HVCs, since they provide "H2-clean" sight lines with minimal contamination. 
In the future, with FUSE or a next-generation satellite with FUV capability, these high-latitude 
sight lines through the H2 holes should be prime targets for studies of HVCs, gas at the disk-halo 
interface, and the IGM. They can be used to address several important astronomical issues: (1) A 
possible correlation of H2 holes with enhanced O VI absorption; (2) The correlation of H2 with H I 
and infrared cirrus structures; (3) High-latitude probes of HVCs, IVCs, and Galactic fountain gas. 

As discussed in § 3.4, the distribution of molecular fractions, /h2, and the level of rotational 
excitation to levels J > 2, can be used to infer the FUV radiation field (H2 pumping rate (5) and 
gas density in the absorbers. Because /h2 oc (n#i?//3) in formation-destruction equilibrium, the 
high-latitude sight lines with high /h2 require either reduced (3 or an enhanced formation rate, tihR- 
The magnitude of the H2 pumping rate ((3) was assessed independently by examining populations 
of the high- J rotational states. Figures 9 and 10 show that the ratios, N(3)/N(l), N(4)/N(0), and 
N(4)/N(2), for high-latitude sight lines with Nh2 > 10 18 cm~ 2 are not greatly different from those 
in the Galactic disk. Thus, it is difficult to argue that the FUV radiation field at z < 200 pc is 
significantly lower than in the disk. Within a factor of two, the absorbers probed in this survey 
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seem to have f3 ~ (3q, the mean value in the Galactic disk. More distant clouds (z > 1 kpc) may 
experience a lower FUV radiation field, as they are farther from OB stars in the Galactic disk. 

Higher H 2 formation rates could arise from more efficient grain catalysis (higher R arising 
from increased grain/gas ratio) or from denser gas (higher n#). At z < 200 pc, the dust /gas 
ratio probably tracks the metallicity and is unlikely to differ much from that in the disk. It is 
more plausible that clouds at the disk-halo interface have been compressed dynamically and have 
densities larger, on average, than quiescent diffuse clouds in the disk. We therefore favor a model in 
which the observed shift in the distribution of fm with Nh arises from more efficient H2 formation 
rate in high-density cirrus clouds. 

In a separate paper (Gillmon & Shull 2005) we identified a correlation between high-latitude 
H2 and infrared cirrus. This correlation means that 100 /tm cirrus maps (Schlegel, Finkbeiner, & 
Davis 1998) can be used to identify the best regions to explore 1ow-Nh2 sight lines in future FUV 
studies. This technique also allows us to characterize the physical properties and total molecular 
mass of the halo clouds. Models of the rotational excitation can be used to estimate the FUV 
radiation field and gas densities, while the IR cirrus can be used to constrain grain temperatures 
(60/100 ^m ratios), define the cloud spatial extent, and characterize the differences from clouds in 
the Galactic disk. By measuring the gas metallicity and dust depletion from UV resonance lines, 
one may even be able to draw connections between the H2, gas metallicity, and dust heating rate. 

We thank James Green, Bart Wakker, Blair Savage, and Ken Sembach for useful discussions, 
and Matthew Browning for his assistance with the modeling of H2 rotational populations. We are 
grateful to the referee, Phillipp Richter, for valuable comments that helped clarify the importance 
on H2 in IVCs and the disk-halo interface. This work was based in part on data obtained for 
the Guaranteed Time Team team by the NASA-CNES-CSA FUSE mission operated by the Johns 
Hopkins University. Financial support to U.S. participants has been provided by NASA contract 
NAS5-32985. The Colorado group also received FUSE support from NASA grant NAG5-10948 for 
studies of interstellar H2. 
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Table 1. GALACTIC HALO TARGET LIST 



Name 


I 

(deg) 


b 
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B 
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Observation ID 
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19040104 
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3 
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2 
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2 
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Table 1 — Continued 



Name I b Type B E(B - V) a Observation ID i oxp b S/N c 
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P1011901 


39.3 


8 














P1011903 


20.8 
















S6011001 


58.5 
















S6011002 


36.9 




PG 0844+349 


188.56 


37.97 


Seyl 


14.83 


0.037 


P1012002 


31.7 


4 


PG 0953+414 


179.79 


51.71 


QSO 


15.29 


0.013 


P1012201 


36.2 


5 














P1012202 


38.1 




PG 1116+215 


223.36 


68.21 


QSO 


14.85 


0.023 


P1013101 


11.0 


6 














P1013102 


10.8 
















P1013103 


8.0 
















P1013104 


10.9 
















P1013105 


33.5 




PG 1211+143 


267.55 


74.32 


Seyl 


14.46 


0.035 


P1072001 


52.3 


5 


PG 1259+593 


120.56 


58.05 


QSO 


15.84 


0.008 


P1080101 


52.3 


7 














P1080102 


48.2 
















P1080103 


51.8 
















P1080104 


106.3 
















P1080105 


104.4 
















P1080106 


62.9 
















P1080107 


95.1 
















P1080108 


32.6 
















P1080109 


29.8 
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Table 1 — Continued 



Name 


I 

(deg) 


b 

(deg) 


Type 


B 


e{b - vy 


Observation ID 


t b 

f-cxp 

(ks) 


S/N c 
(pixel- 1 ) 


PG 1302-102 


308.59 


52.16 


QSO 


15.18 


0.043 


P1080201 


31.9 


4 














P1080202 


32.3 
















P1080203 


81.9 




PKS 0405-12 


204.93 


-41.76 


QSO 


15.09 


0.058 


B0870101 


(71.1) 


4 


PKS 0558-504 


257.96 


-28.57 


QSO 


15.18 


0.044 


P1011504 


45.1 


5 














C1490601 


(48.0) 




PKS 2005-489 


350.37 


-32.60 


BLLac 


13.40 


0.056 


P1073801 


12.3 


5 














C1490301 


(22.5) 
















C 1490302 


(13.3) 




PKS 2155-304 


17.73 


-52.25 


BLLac 


13.36 


0.022 


P1080701 


19.2 


10 














P 1080703 


65.4 
















P 1080705 


38.6 




Ton S180 


139.00 


-85.07 


Seyl.2 


14.60 


0.014 


P1010502 


16.9 


4 


Ton S210 


224.97 


-83.16 


QSO 


15.38 


0.017 


P1070302 


41.0 


5 


VII Zw 118 


151.36 


25.99 


Seyl 


15.29 


0.038 


P1011604 


67.2 


4 



P1011605 25.9 
S6011301 49.9 



& E(B — V) values are from Schlegcl et al. (1998) by way of NED (http://nedwww.ipac.caltech.edu). Since 
these values are inferred from IRAS dust maps for a nonredshifted elliptical galaxy, they may not be appro- 
priate for other source types. Thus, these values should be treated as an indication of relative dust column 
along the sight line as opposed to the actual reddening for these AGN sources. 

b Exposure times are estimates of the actual post-calibration exposure times. The exposure times were 
partially corrected by subtracting the duration of the subexposures during which the detector was off. Because 
the calibration may omit additional subexposures with low flux or burst activity, the listed exposure times 
are upper limits on the actual exposure times. For those targets without available subexposure information, 
the full uncorrected exposure time is given in parentheses. 

c For targets with more than one observation, the S/N for the coadded dataset (including all listed observa- 
tions) is given in the row with the target name. 
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Table 2. ROTATIONAL LEVEL POPULATIONS 



Name 


log N(0) 


log N(l) 


log N(2) 


log N(3) 


log N(4) 


log N(5) 


^dopp 




(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 


(cm" 2 ) 


(km s- 1 ) 



3C 249.1 


18.40 ±8;i 


18-84 ±g : }4 


17 


.03 


±8: 


80 
98 


16 


41 


±J: 


19 
55 








- 


6.1 ± 2 :i 


3C 273 (v = -5) 


<13.51 


14.30 ±8:81 
























- 


linear 


(v = 26) 


14.89 ±g;» 


15.50 ±g fi 


14. 


,87 


±1: 


10 

.07 


14 


.71 


±8: 


07 
05 








- 


4 Q +0.9 


ESO 141-G55 


18.90 ±8;}} 


i9.io ±8:81 


17 


.34 


±8 


44 
.79 


16 


10 


±,V 


26 
52 


14.64 


± 


1.56 
0.25 


- 


3.8 ±\\ 


H 1821+643 


16.97 ±0; 2 6 


i7.7i ±8:81 


17. 


.18 


±8. 


Ifi 
2?, 


16 


,78 


+° 
± o 


31 
61 








- 


i.7 ±8:? 


HE 0226-4110 b 


<13.90 


<14.06 
























- 


- 


HE 1143-1810 


15.91 ±J;|g 


16.37 ±J:18 


15. 


.35 


±0 


34 
31 


14 


,93 


±8: 


97 
25 








- 


7.0 ±§;6 


HS 0624+6907 


19.39 ±°H 


19.60 ±8:11 


17. 


.95 


±?: 


.31 
56 


17 


.13 


±?: 


57 
13 


15.49 


± 


0.33 
0.37 


14.63 ±8:81 


e.s ±f : 2 


MRC 2251-178 


<14.12 


14.54 ±8' 2 ? 
























- 


linear 


Mrk 9 


18.87 ±g;18 


19.18 ±8i} 


17. 


.13 


±° 


.35 
46 


15. 


.40 


±° 
^0. 


20 
20 








- 


5.4 ±8:| 


Mrk 106 


15.95 ±8;i 


15.82 ±8: 2 | 


14 


.98 


+ (l 


14 
11 


14 


.73 


+° 


16 
15 








- 


14.9 ±1:1 


Mrk 116 


18.82 ±8:11 


18.73 ±8'?? 


17 


.12 


±1! 


47 

s;> 


16 


.89 


±1! 


58 
71 


14.76 


± 


0.40 
0.24 


- 


4-5 ±U 


Mrk 205 


15.63 ±8;«8 


16.39 ±8ig 


15 


.44 


±8 


00 

.28 


15 


.40 


±8: 


00 
38 








- 


7-0 ±5:8 


Mrk 209 


<14.10 


<14.25 
























- 


- 


Mrk 290 


15.80 ±J : §« 


15.78 ±§ : §| 


15. 


26 


+ 2 


21 
27 


15 


03 


+ 2 ' 


28 
26 








- 


8-7 ±1:5 


Mrk 335 


18.44 ±0-11 


18.59 ±8:8? 


16. 


,78 


+° 


50 
44 


16 


.07 


± o. 


51 

24 


14.23 


± 


0.15 
0.12 


- 


s.s ±?:l 


Mrk 421 


13.89 ±8;i 


14 40 +0- 13 


14 


.01 


+ (l 
± () 


25 
25 
















- 


linear 


Mrk 478 


<14.18 


<14.33 
























- 


- 


Mrk 501 


<14.20 


14-78 ±8 : }g 
























- 


linear 


Mrk 509 


17-34 ±8;1| 


17.69 ±8: 2 ^ 


16 


.37 


+° 
± () 


62 
fil 


15 


.63 


_i_0. 

± () 


38 
26 








- 


5.8 ±i:§ 


Mrk 817 


<13.64 


<13.80 
























- 


- 


Mrk 876 (v = -33) 


15.18 ±8:11 


15.53 ±g: 22 


14. 


61 


+ 2 
± () 


22 
2fi 


14 


50 


+ 2 ' 
± (). 


15 

25 










5.7 ±iv 


(v = -3) 


15.81 ±lf 


16.41 ±J;98 


15 


.59 


+ 2 


11 

38 


15 


.22 


+ 2 - 


23 
30 








- 


7-3 ±1:1 


Mrk 1095 


18.24 ±8;|| 


18.58 ±8:1 


17. 


.10 


± () 


77 
72 


16 


.49 


± () 


92 
50 


14.90 


± 


0.22 
0.18 


- 


9.2 ± 2 : ! 


Mrk 1383 


<13.77 


14.35 ±S:ig 
























- 


linear 


Mrk 1513 


16-12 ±hii 


16.02 ±i : 0| 


15. 


.23 


±8 


54 
23 


14 


,87 


±8: 


27 
19 








- 


11-7 ±3-5 


MS 0700.7+6338 


i8.4i ±8:18 


18.46 ±8:|| 


17 


.27 


±? 


86 

.57 


15. 


37 


±J 


75 
41 










7-6 ±1:1 


NGC 985 


15.57 ±Hi 


15.80 ±tf 8 


14. 


93 


±5 


86 
2(1 


14 


,67 


±,V 


06 
20 










8.1 ±1:8 


NGC 1068 


17.84 ±8;8| 


17.82 ±8:?f 


15 


.84 


+° 


66 
40 


15 


.39 


+° 


60 
26 










3 4 +°- 8 

^o.s 


NGC 1705 


<13.78 


<13.94 




























NGC 4151 


15.95 ±8;I 2 


16.55 ±8:11 


15 


.57 


+° 


.57 
12 


15 


19 


_|_0, 


26 
11 










7 2 + 1 - 1 


NGC 4670 


14.31 ±8|t 


14.51 ±8:11 


























linear 


NGC 7469 


19-41 ±8:58 


19.32 ±8:81 


17. 


,77 


+° 


24 
.35 


16 


39 


+° 

^0 


75 
35 


14.90 


± 


0.27 
0.14 




5.5 ±J'5 


PC 0804+761 


18.08 ±8;Jg 


18.52 ±8'?f 


16 


63 


+ (l 

^0 


54 
37 


15 


.86 


^0 


32 
16 










8 .7±i; 2 


PG 0844+349 


17.64 ±8;i 


18.09 ±° o: \l 


16. 


,04 


+° 


.78 
40 


15 


.04 


+° 
± () 


26 
21 










4-7 ±i;g 


PG 0953+414 


14.10 ±8;i 


14.76 ±8: J? 


14 


09 


+ (l 
± () 


25 
25 


14 


.39 


+° 
± () 


23 
20 










linear 


PG 1116+215 


<13.78 


<13.93 




























PG 1211+143 


17.80 ±8;Ji 


18.23 ±8:81 


16. 


,72 


+° 


57 
59 


15 


.77 


+° 
± () 


81 

25 










4.3 ±1:8 


PG 1259+593 


13.77 ±8; 2 1 


14.34 ±8:g 


14 


,06 


+° 


25 
25 


14. 


.23 


_|_0, 


25 
25 










linear 


PG 1302-102 


14.79 ±8:1? 


15.37 zLq 23 


14. 


,84 


+ (l 
± () 


50 
19 


14 


.71 


+° 
± () 


43 
18 










11.0 ±1:8 


PKS 0405-12 


14.91 ±8;g 


15.56 ±8 : ?I 


14. 


93 


+° 
± () 


10 
08 


14 


.96 


+° 
± (). 


13 
10 










11.8 ±1:5 


PKS 0558-504 


i4.5i ±8;3:2 


15.22 ±8:11 


14 


,65 


+° 


12 
11 


14. 


.52 


+° 


16 
16 










10.5 ±|:| 


PKS 2005-489 


14.32 ±8|1 


14.74 ±8:5| 


14 


,36 


+ (l 
± () 


25 
25 


14 


.29 


+° 
± () 


25 
25 










linear 


PKS 2155-304 


13.59 ± : 25 


14.04 ±8:il 


























linear 


Ton S180 


<13.98 


<14.14 
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Table 2 — Continued 



Name 


log N(0) 
(cm -2 ) 


log N(l) 
(cm -2 ) 


log N(2) 
(cm -2 ) 


log N(3) 
(cm -2 ) 


log N(4) 
(cm- 2 ) 


log N(5) 
(cm- 2 ) 


^dopp 

(km s— x ) 


Ton S210 
VII Zw 118 


16.00 ±tfl 
18.38 ±q \1 


16.44 ±l;0i 
18.65 ±° ° l 


16.22 ±°f 4 


16.01 ±g;Si 


14.40 ±8;if 




4.4 ±H 

6.5 ±i;? 



a Entry of "linear" means that lines fall on the linear part of the curve of growth, and thus there is no constraint 
on dopplcr b. 

b Savagc ct al. (2005) have detected H 2 in this sightline, with log N(0) = 13.85 ± 0.3, log N(l) = 14.06 ± 0.3, 
log N(2) = 13.44±°J, and log N(3) = 13.5lt?;|. 
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Table 3. H2 and H I Velocity Components (LSR) 



Name 


V[ri L) 




Umax { tl I) 


V (tl2 ) 


IVLj, rl2 status 






(Km s ) 


^Kin s ) 


(Km s ) 


firm 1*1 

[Km s ) 


[Km s ) 


in otes 


3C 249.1 


-50, -20, (1,9) 


9 


-50, -20,1,9 


12 


—50, Possibly 


2 


3C 273 


(-13, -6), 0, 25 


-6 


-13, -6, 


-5 


25, Yes 


1 




-13,-6,0, (25) 


25 


25 


26 


25,Yes 


3 


ESO 141-G55 


-45,-23, (-1,2), 30 


-1 


-45, -23, -1, 2, 30 


-3 


—45, None 


2 


H 1821+643 


-87, -23, (-10. -3) 


-3 


-23, -10, -3 


-3 






HE 0226-4110 














HE 1143-1810 


(-6, -5,4), 13,56 


4 


-6, -5,4, 13 


3 






HS 0624+6907 


-27, (-10, -3,0, 10) 


10 


-27,-10,-3,0, 10 


6 






MRC 2251-178 


-8, (0), 59 


-8 


-8,0 


-8 






Mrk 9 


-40, (-8, -3,-1,4) 


-8 


-40,-8,-3,-1,4 


-2 


—40, Possibly 


2 


Mrk 106 


-40, (-6, -2, -2) 

7 \ 1 J / 


-6 


-6, -2, -2 


-6 


-40, Possibly 


1 


Mrk 116 


-39, -16, (1, 7) 


7 


-39, -16, 1, 7 


4 


-39, Possibly 


2 


Mrk 205 


-48, (-4, 2) 


2 


-4,2 


3 


-48, Possibly 


1 


Mrk 209 










—55, None 














-100, None 




Mrk 290 


f — 6, — 5 , — 5f>) s 

y ") ■-'lit u 1 


-5n 


-6 -5„ -5 6 


-8 






Mrk 335 


-27, (-9, -3, -2) 


-2 


-27, -9, -3, -2 


-1 


-27,Possibly 


2 


Mrk 421 


-61, (-23, -11,-8) 


-23 


-23,-11,-8 


-15 


-61, None 


1 


Mrk 478 














Mrk 501 


-83, -38, (-3, -1,2) 


2 


-3, -1,2 


-1 


-38,None 


1 


Mrk 509 


(-7, 6), 60 


-7 


-7,6 


2 


60,Yes 


1 


Mrk 817 










-40, None 




Mrk 876 


-52, (-30), -7 


-30 


-30 


-33 


-30,Yes 


3 




-52, -30, (-7) 


-7 


-7 


-3 


-30,Yes 


1 


Mrk 1095 


-13, (4,9) 


9 


-13,4,9 


16 






Mrk 1383 


(-7,-2) 


-2 


-7, -2 









Mrk 1513 


-29, (1, 5, 10) 


10 


-29, 1,5, 10 


6 


-29, None 


2 


MS 0700.7+6338 


-11, (-6, -1,6) 


-6 


-11, -6, -1, 6 


6 






NGC 985 


(-8,2) 


2 


-8,2 


-7 


_ 


_ 


NGC 1068 


-58, -6, (5) 


5 


-6,5 


16 


-58, None 


1 


NGC 1705 










87,Nonc 




NGC 4151 


(-41, -29), -24, -21,0 


-41 


-41,-29,-24,-21 


-35 


-29, Yes 


3 


NGC 4670 


-28, (-19, -8) 


-8 


-19, -8 


-11 






NGC 7469 


-13, (-3,0) 





-13, -3,0 


-1 






PG 0804+761 


-58, (-10, -5,6) 


-5 


-58, -10, -5,6 


4 


-58,Yes 


2 


PG 0844+349 


-12, (1,6) 


6 


1,6 


8 






PG 0953+414 


-49, (-4) 


-4 


-4 


3 


-49, None 


1 


PG 1116+215 










-42,Yes 




PG 1211+143 


-35, (-20, -9) 


-20 


-35,-20,-9 


-12 


—35, None 


2 


PG 1259+593 


-54, (-6) 


-6 


-6 


-5 


-54,Yes 


1 


PG 1302-102 


-53, (-22, -3, -1), 14 


-1 


-22,-3,-1,14 


-8 






PKS 0405-12 


(-13, 1,4), 18 


1 


-13, 1,4, 18 









PKS 0558-504 


(4,6),26,27,71 


6 


4,6,26,27 


10 


26,None 


2 












71, None 


1 


PKS 2005-489 


-26, (2) 


2 


2 


-1 
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Table 3 — Continued 



Name 


v(H I) a 


fmin(H I) b 


»max(H I) c 


f(H 2 ) d IVC, H 2 status 


IVC f 




(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) (km s- 1 ) 


Notes 


PKS 2155-304 


(-4), 14 


-4 


-4 


-8 




Ton S180 












Ton S210 


-30, (-10, -6) 


-6 


-30,-10,-6 


-10 




VII Zw 118 


(-7, -3,7) 


7 


-7,-3,7 


-2 





a «(H I) lists the 21 cm velocity components fitted by Wakker et al. (2003), excluding high-velocity 
components at \v\ > 90 km s _1 . The components in parentheses are those we have selected as most likely 
to correspond to the H2 absorption, based on the coincidence of radial velocities (§ 2.3). The sum of the 
column densities for the components in parentheses is the value we adopt for Njji- 

b u m j n (H I) indicates the velocity component with the mimimum column density that, based on the coin- 
cidence of radial velocities, could be solely associated with the H2. The column density of this component 
provides the lower bound for Njji. 

c fmax(H I) lists all velocity components that, based on the coincidence of radial velocities, could feasibly 
be associated with the H2 . The sum of the column densities of these components provides the upper bound 
for Nhi- 

d In general, the uncertainty on f(H2) is ~ ±5 km s —1 for (S/N)t,i n > 4. 

c Velocity and H2 status of IVCs presented in Table 2 of Richter et al. (2003). 

'Notes on IVCs: (1) IVC does not contribute to fitted H2 lines; (2) Possible IVC contribution on edge 
of fitted H2 lines; (3) IVC is main fitted H2 component. Although we generally did not measure H2 in 
IVCs, we did so for targets noted as (3). For 3C 273, we do not consider component at 25 km s _1 to be 
true IVC. For Mrk 876, we fitted the IVC lines blended with LSR lines (double component fits). For NGC 
4151, the H2 lines may have some contributions from non-IVC gas. 

g Since this sight line has two H I components at the same velocity, —5 km s _1 , we distinguish them by 
"b" for broad and "n" for narrow. 
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Table 4. H 2 SIGHTLINE SUMMARY 



Name 


log N H2 


log N H i 


beam a 


log f H 2 


T01 


T'exc 




(cm- 2 ) 


(cm- 2 ) 


size 






(K) 


(K) 


3C 249.1 


18.98 ±g;i| 


i n on 

20.25 ±8;|2 


3 




-1.01 


144 ±Jf 4 


406 ±317 


3C 273 (v = -5) 


14 -30 ±8:81 


20.14 ±8;§ 7 


3 




-5.54 






(v = 26) 


15.72 ±0; 2 | 


19.42 ±8;8| 


3 




-3.40 


216 ±^ 24 


645 ±71 


ESO 141— G55 


-i o oo i 07 

19.32 ±8;8f 


OO >~70 1 08 

20.70 ±8 67 


2 




1.12 


98 ±H 


416 ±72 


H 1821+643 


1 — o-i I 1 1 

17.91 ±g;i3 


20.43 ±8;^ 


4 




-2.22 




494 ±168 


HE 0226— 4110 b 


<14.29 


>19.50 


2 


< 


-4.92 






HE 1143—1810 


-i r - ^ I 1 1Q 

16 - 54 ±0:68 


oo a '-j i n oi 
20.47 ±i 19 


3 




3.63 


-1 rrn 1 206 
150 ±74 


484 ±185 


Hb 0624+6907 


-1 n oo 1 1 

19.82 ±8;i8 


o o o o i 9 

20.80 ±i 06 


4 




0.76 


100 ±f| 


498 ±34 


MHO 2251 — 178 


i/i k /i _l0 28 

14.54 ±q'i7 


oo no 1 01 

20.39 ±\»l 


3 




-5.55 






Mrk 9 


i n o/^ i 09 

19.36 ±8;8l 


oo n a 1 04- 

20.64 ±i "| 


3 




1.02 


H5 ±25 


215 ±37 


MrR luo 


1 R oo _i_0.21 
16.23 ± 15 


on oc 4-0.03 
zu.oo zn 83 


A 
1 




-3.82 


Ra 4-21 
68 ± 12 


C7Q 4_ 1 OQ 

oto ±lzo 


Mrk 116 


1 n no i fl 11 

19.08 ±8;ig 


oo a 1 i n 07 
20.41 ±8;^ 


4 




-1.06 


t-i 1 1 fi 

71 ±J6 


437 ±82 


Mrk 205 


16.53 ±°;13 


oo .i o i n oi 

20.40 ±8;8| 


4 




-3.57 




762 ±542 


MrR zuy 


\14.4o 


^ iy. / v> 


Q 
O 


< 


-4.95 






Mrk 290 


16.18 ± 2 : °i 


oo -it i n oi 
20-11 ±0:93 


4 




-3.62 


76 ±|| 


592 ±260 


T\T„l r not; 

Mrk ooo 


1 a _i_0.08 
18.83 ± . 08 


on /i q 1 0.14 
20.43 ± 1 30 


3 




1.32 


oo J_27 
92 ± 14 


A Q A 4_ £R 

434 ±56 


Mrk 421 


14.63 ±8;^ 


m nj 1 03 

19.94 ±1 16 


3 
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Table 4 — Continued 



Name 


log N H 2 


log N H i 


bcam a 


log fH2 


Toi 


^exc 




(cm -2 ) 


(cm -2 ) 


size 




(K) 


(K) 


Ton S210 


16.57 ±i;i| 


20.19 ±°/°j 


4 


-3.32 


144 ±l| 2 




VII Zw 118 


18.84 ±°;}° 


20.56 i&gg 


3 


-1.43 


108 ±|| 


552 ±80 



a Source of 21 cm data: (1) Leiden-Dwingcloo Survey (35' beam), (2) Villa Elisa tele- 
scope (34' beam), (3) Green Bank 140-ft telescope (21' beam), (4) Effelsberg telescope 
(9'. 7 beam). 

b Savagc et al. (2005) have detected H 2 in this sightline, with log N H2 fts 14.39 ±0.43. 



- 28 - 



NGOKTO 




Fig. 1. — Aitoff projection showing locations of 45 AGN targets from our high-latitude H2 survey. 
In the northern Galactic hemisphere, we identified 8 high-latitude (b > 54°) sight lines with low 
Nh2 < 10 15 ' 1 cm -2 , color-coded in red. These targets suggest a widespread H2 hole (60° < I < 180° 
and b > 54°). Sightlines with higher Nh2 are labeled in green and blue. We sampled targets with 
Galactic latitude \b\ > 20° to minimize absorption by interstellar gas in the Galactic disk. The 
low column density "H2 Holes" at |6| > 54° may be related to the "Northern Chimney" (region of 
low Na I absorption, Lallement et al. 2003) and the "Lockman Hole" (region near of low Nhi near 
i = 135° - 175° and b = 46° - 60°, Lockman et al. 1986, see also Figure 10). These "H 2 -clean" 
sight lines are ideal for FUV studies of the IGM and HVCs. 
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Fig. 2. — Top Panel: 21 cm spectrum of Mrk 509 from Leiden-Dwingeloo Survey. Middle and 
bottom panels: FUSE spectra binned by 4 pixels. Dotted lines indicate the three H I components 
identified by Wakker et al. (2003) listed in top panel. FUSE data have been shifted to align the H I 
emission and Ar I absorption. A comparison of H I and H2 shows that either the — 7 km s -1 H I 
component or the 6 km s _1 H I component, or both, could be associated with the H2 absorption. 
We adopt the sum of the two column densities, log Nh2 = 20.58, for Nhi- Since either component 
could be solely responsible for the H2 absorption, we choose the lower column density, log Nh2 
=20.00, for our lower limit on Nhi- Richter et al. (2003) detect H2 in the 60 km s _1 component, 
but this absorption does not contribute to our fitted H2 lines. Since no other H I components can 
contribute to the H2 absorption, we adopt the sum of the —7 km s _1 component and 6 km s _1 
component (plus a systematic uncertainty, A log Nhi = 0.03) for the upper limit on Nhi- 
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Fig. 3. — Upper two panels: FUSE spectrum of PG 1211+143, showing expected locations of (0-0) 
through (17-0) H2 Lyman bands (red) and (0-0) through (4-0) Werner bands (green). This sight 
line has log Nh2 = 18.38lo' 14 in J = — 3. Lower panel: Blowup of the (4-0) Lyman band of H2, 
showing R(0), R(l), P(l), R(2), P(2), and R(3) absorption lines from J = 0, 1, 2, and 3 together 
with interstellar lines of Ar I A1048.220 and Fe II A1055.262. 
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Fig. 4. — Examples of FUSE spectra of four AGN, showing increasing H2 column densities from 
the Lyman (4-0) band. First panel: Ton S180 shows no detectable to a limit log Nh2 < 14.37 
(< 13.98 in J = and < 14.14 in J = 1). Second panel: PKS 2005-489 with log N H2 = 15.07+ ' 10 



Third panel: NGC 4151 with log N H2 = 16.70+^. Fourth panel: ESO 141-G55 with log N H2 
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Fig. 5. — Distribution of Nh2 (the LSR absorption only) for 38 of our 45 AGN sight lines. An 
additional 7 sight lines (16% of the sample) showed no detectable low- velocity H2 above a threshold, 
log N H2 = 13.8-14.6, depending on the S/N (2-11 per pixel) and resolution (R = 15,000 - 20,000 
for FUSE). The dip in numbers between log Nh2 = 17-18 may divide two populations (see also 
Fig. 9), although this column density bin is sensitive to CoG uncertainties and is based on small 
numbers of sight lines. 
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Fig. 6. — Transitions in molecular fraction, /h2 5 with total hydrogen column density, Nh = (Nhi 
+ 2Nh2)- In our FUSE survey of 45 high-latitude sight lines, the transition occurs at log Nh ~ 
20.38 ±0.13, approximately a factor of two below the transition (log Nh = 20.7) in the Copernicus 
survey (Savage et al. 1977). This shift to lower Nh could arise from reduced photodissociation 
(lower UV field), from enhanced H2 formation rates (higher tlh in compressed cirrus clouds) or 
both. 



o 



500 



400 



300 



200 



100 



MEAN 



L_ 
19.0 



□ 



c] 



■ B 



i: 



□ 



□ 



□ [] 



logN(H 2 ) > 17 



[] 



[] 



19.5 20.0 20.5 

log[N(HI) + 2N(H 2 )] cm" 2 



21.0 



Fig. 7. — Mean rotational temperature in our Galactic high-latitude survey is (Tqi) = 124 ±8 K 
(median 121 K). For 15 sight lines with log Nh2 > 17 (solid squares), we find (Tqi) = 109 ± 7 K. 
These temperatures are somewhat higher than those found in the Copernicus H2 survey (77 ± 17 K, 
Savage et al. 1977) and in our FUSE Galactic disk survey (86 ± 20 K, Shull et al. 2005). The 
range of Tqi is 68-252 K, consistently higher than the Galactic disk values, all but two of which lie 
between 55-120 K. 
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Fig. 8. — The mean excitation temperature for J > 2 is (T cxc ) = 505 ± 28 K (median 501 K), 
based on 30 sight lines with detectable J = 3 lines (in some cases J = 4 and J = 5). This value is 
slightly larger than the mean, 326 ± 125 K, seen in the FUSE disk H2 survey (Shull et al. 2005). 
The distribution of T cxc is a measure of non-thermal excitation by UV pumping and perhaps H2 
formation on grain surfaces (Browning et al. 2003). 



-36- 



22 ; 

20 - 



I 18 



D) 



J3 _ • Halo Survey TP 
ID Disk Survey » 

14 r , , , , . . . 

-5 -4 -3 -2 -1 

logN(3)/N(1) 





Fig. 9. — Ratios of rotational-state column densities for J = 3/1, J = 4/2, J = 4/0 vs. Nh2 for 
sight lines in our surveys in the Galactic halo (blue) and Galactic disk (green). Of the 45 halo sight 
lines, 8 have measured column densities for both J = 3 and J = 4, while 32 have J = 3 but lack 
detectable J = 4 because of insufficient total Nh2- For all sight lines with log Nh2 > 18, the degree 
of rotational excitation is comparable (Shull et al. 2005), except for a population of 16 halo sight 
lines with low Nh2 and highly excited J = 3 [log N(3)/N(l) > —1.5]. A comparison of the high-Nn2 
data with models of H2 formation and rotational excitation (Browning et al. 2003) suggests that 
the UV radiation fields in the low Galactic halo and disk are similar. To produce the observed shift 
to lower Nh in the transition of /h2; the H2 formation rate may be enhanced in the halo absorbers, 
because of higher uh in compressed cirrus clouds. 
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Fig. 10. — Models of H2 rotational excitation ratios, N(J), using the code of Browning, Tumlinson, 
& Shull (2003). Plots show ratios, JV(4)/JV(2) and N(3)/N(l), versus N H2 , for data (blue points, 
same as in Fig. 9) and for a grid of models (green points). Models assume the mean Galactic H2 
formation rate (R = 3 x 10 -17 cm 3 s _1 ), and a range of parameters, FUV radiation fields, varying 
by a factor of two above and below the mean ISM value (Jfuv = 2 x 10~ 8 photons cm' 2 s' 1 
Hz -1 between 912-1120 A). The model grid ponts include variations in cloud temperatures (T = 
20-150 K) and density (nn = 5-800 cm' 3 ). Note the good agreement between data and models, 
which show increasing excitation, N(3)/N(l), for "thinner" clouds, log Nh2 < 18. 
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Fig. 11. — Locations in northern Galactic hemisphere (b > 45°) of 14 high-latitude AGN sight lines 
(white circles). The symbols inside white circles indicate log Nh2 in the ranges: 14.0-15.1 (Star); 
15.1-17.0 (Diamond); 17.0-20.0 (Square). The targets sample the distribution of Nhi in the Leiden 
Dwingeloo Survey (Hartmann & Burton 1997), colored as deep-blue (log Nhi = 19.7-20.0) or light- 
blue (20.0-20.1). Gas at higher Nhi is shown in green, yellow, or red. The extended Lockman Hole 
(log N H i « 19.5-20.0) is seen in blue/purple {£ = 135° - 175° and b = 46° - 60°) at lower right. 



